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Abstract 


In his Habilitationsvortrag, Riemann described infinite dimensional manifolds pa¬ 
rameterizing functions and shapes of solids. This is taken as an excuse to describe 
convenient calculus in infinite dimensions which allows for short and transparent 
proofs of the main facts of the theory of manifolds of smooth mappings. Smooth 
manifolds of immersions, diffeomorphisms, and shapes, and weak Riemannian met¬ 
rics on them are treated, culminating in the surprising fact, that geodesic distance 
can vanish completely for them. 
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1 Introduction 

At the very birthplace of the notion of manifolds, in the Habilitationsschrift EO] 
end of section I], Riemann mentioned infinite dimensional manifolds explicitly: 

“Es giebt indess auch Mannigfaltigkeiten, in welchen die Ortsbestim- 
mung nicht eine endliche Zahl, sondern entweder eine unendliche Reihe 
oder eine stetige Mannigfaltigkeit von Grossenbestimmungen erfordert. 
Solche Mannigfaltigkeiten bilden z.B. die moglichen Bestimmungen 
einer Function fiir ein gegebenes Gebiet, die moglichen Gestalten einer 
raumlichen Figur u.s.w.” 

The translation into English in reads as follows: 

“There are manifoldnesses in which the determination of position re¬ 
quires not a finite number, but either an endless series or a continuous 
manifoldness of determinations of quantity. Such manifoldnesses are, 
for example, the possible determinations of a function for a given re¬ 
gion, the possible shapes of a solid figure, &c.” 

If one reads this with a lot of good will one can interpret it as follows: When 
Riemann sketched the general notion of a manifold, he also forsaw the notion of an 
infinite dimensional manifold of mappings between manifolds, and of a manifold of 
shapes. He then went on to describe the notion of Riemannian metric and to talk 
about curvature. I will take this as an excuse to describe the theory of manifolds 
of mappings, of diffeomorphisms, and of shapes, and of some striking results about 
weak Riemannian geometry on these spaces. See [7] for an overview article which 
is much more comprehensive for the aspect of shape spaces. 

An explicit construction of manifolds of smooth mappings modeled on Frechet 
spaces was described by [TD]. Differential calculus beyond the realm of Banach 
spaces has some inherent difficulties even in its definition; see section [5J Smooth¬ 
ness of composition and inversion was first treated on the group of all smooth 
diffeomorphisms of a compact manifold in |41) : however, there was a gap in the 
proof, which was first filled by m- Manifolds of C^-mappings and/or mappings 
of Sobolev classes were treated by [35], [13], Smale-Abraham [T], and [S3]. Since 
these are modeled on Banach spaces, they allow solution methods for equations 
and have found a lot of applications. See in particular m- 


2 A short review of convenient calculus in infinite 
dimensions 

Traditional differential calculus works well for finite dimensional vector spaces and 
for Banach spaces. Beyond Banach spaces, the main difficulty is that composition 
of linear mappings stops to be jointly continuous at the level of Banach spaces, for 
any compatible topology. Namely, if for a locally convex vector space E and its 
dual E' the evaluation mapping ev : A x A' —>■ R is jointly continuous, then there 
are open neighborhoods of zero U C E and U' C E' with ey{U x U') C [—1,1]. 
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But then U' is contained in the polar of an open set, and thus is bounded. So E' 
is normable, and a fortiori E is normable. 

For locally convex spaces which are more general than Banach spaces, we sketch 
here the convenient approach as explained in and [33] . 

The name convenient calculus mimicks the paper |()3j whose results (but not 
the name ‘convenient’) was predated by [12], m. [I4] . They discussed compactly 
generated spaces as a cartesian closed category for algebraic topology. Historical 
remarks on only those developments of calculus beyond Banach spaces that led to 
convenient calculus are given in [531 end of chapter I, p. 73ff]. 

2.1 The c°°-topology 

Let E he a locally convex vector space. A curve c : K ^ A is called smooth 
or C°° if all derivatives exist and are continuous. Let C°°{M.,E) be the space of 
smooth curves. It can be shown that the set C°°(]R, A) does not entirely depend 
on the locally convex topology of E, only on its associated bornology (system of 
bounded sets); see [351 2.11]. The final topologies with respect to the following 
sets of mappings into E coincide; see [331 2.13]: 

1 . C°^{R,E). 

2. The set of all Lipschitz curves (so that : t ^ s, |t|,|s| < C} is 

bounded in E, for each C). 

3. The set of injections Eb —t E where B runs through all bounded absolutely 
convex subsets in E, and where Eb is the linear span of B equipped with 
the Minkowski functional ||a:||B := inf {A > 0 : a: € XB}. 

4. The set of all Mackey-convergent sequences —)• a; (there exists a sequence 

0 < Xn oo with Xn(xn — x) bounded). 

This topology is called the c°°-topology on E and we write c^E for the resulting 
topological space. 

In general (on the space V of test functions for example) it is finer than the 
given locally convex topology, it is not a vector space topology, since addition is 
no longer jointly continuous. Namely, even c°°{'D xV) ^ c°°'D x c°°T>. 

The finest among all locally convex topologies on E which are coarser than 
c^E is the bornologification of the given locally convex topology. If A is a Frechet 
space, then c°°E = E. 

2.2 Convenient vector spaces 

A locally convex vector space E is said to be a convenient vector space if one of 
the following equivalent conditions holds (called c°°-completeness); see [331 2.14]: 

1. For any c G C"^(R, E) the (Riemann-) integral c{t)dt exists in E. 

2. Any Lipschitz curve in E is locally Riemann integrable. 
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3. A curve c : M —>■ A is C°° if and only if A o c is C°° for all X € E*, where E* 
is the dual of all continuous linear functionals on E. 

• Equivalently, for all A € E', the dual of all bounded linear functionals. 

• Equivalently, for all A € V, where V is a subset of E' which recognizes 
bounded subsets in E. 

We call this scalarwise C°°. 

4. Any Mackey-Cauchy-sequence (i.e., tnm{xn — Xm) 0 for some tnm —>■ oo 
in R) converges in E. This is visibly a mild completeness requirement. 

5. If B is bounded closed absolutely convex, then Eb is a Banach space. 

6. If / : R —> A is scalarwise Lip^, then / is Lip^, for fc > I. 

7. If / : R —)• A is scalarwise C°° then / is differentiable at 0. 

Here a mapping / : R —>■ A is called Lip^ if all derivatives up to order k exist 
and are Lipschitz, locally on R. That / is scalarwise C°° means A o / is C°° for 
all continuous (equiv., bounded) linear functionals on E. 

2.3 Smooth mappings 

Let E, and E be convenient vector spaces, and let U C E be c°“-open. A mapping 
/ : 1/ —>■ E is called smooth or C°°, if / o c S C°“(R, F) for all c G C°°(R, U). See 

[331 3.II]. 

If E is a Frechet space, then this notion coincides with all other reasonable 
notions of C^-mappings; see below. Beyond Frechet spaces, as a rule, there are 
more smooth mappings in the convenient setting than in other settings, e.g., C^. 

2.4 Main properties of smooth calculus 

1. For maps on Frechet spaces this coincides with all other reasonable defini¬ 
tions. On R^ this is non-trivial; see m or [31 3.4]. 

2. Multilinear mappings are smooth iff they are bounded; see [31 5.5]. 

3. It E D U ^ > E is smooth then the derivative df : U x E ^ F is smooth, 
and also df : U ^ L{E, F) is smooth where L{E, F) denotes the space of 
all bounded linear mappings with the topology of uniform convergence on 
bounded subsets; see [31 3.18]. 

4. The chain rule holds; see [31 3.18]. 

5. The space C°°{U,F) is again a convenient vector space where the structure 
is given by the injection 

H C°°(R,R), 

cGC“(R,f7),^GF* 


C°°{U,F) 
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and where M) carries the topology of compact convergence in each 

derivative separately; see [33l 3.11 and 3.7]. 

6. The exponential law holds; see [33j 3.12].: For c°°-open V C F, 

C^{U,C°°{V,G)) ^ C^{U X V,G) 

is a linear diffeomorphism of convenient vector spaces. Note that this (for 
U = 'M.) is the main assumption of variational calculus. Here it is a theorem. 

7. A linear mapping f : E ^ G°°{V,G) is smooth (by (2) equivalent to 

bounded) if and only if E ^ ) G°°iy,G) —y G is smooth for each 

V & V. {Smooth uniform boundedness theorem; see [331 theorem 5.26].) 

7. A mapping f : U ^ L{F, G) is smooth if and only if UL{F, G) > G is 
smooth for each v € F, because then it is scalarwise smooth by the classical 
uniform boundedness theorem. 

8. The following canonical mappings are smooth. This follows from the expo¬ 
nential law by simple categorical reasoning; see [33l 3.13]. 

ev : G°°{E, F) x E ^ F, ev(/, x) = f{x) 
ins : A —>• G°°{F,E x F), iiis{x){y) = {x,y) 

( )'^ :G°^{E,G°°{F,G)) ^ G°°{E X F,G) 

( y ■.G°^{ExF,G)^G°°{E,G°°{F,G)) 
comp : G^{F,G) x G^{E,F) G°°{E,G) 

G°°{ , ) : G°°{F, Fi) X C“(Ai, A) ^ 

^ C“(C“(A,F),C“(Ai,Fi)) 

{f, 9 )^ fo hog) 

n :Y[g^{e,,f,) ^ 


This ends our review of the standard results of convenient calculus. Just for 
the curious reader and to give a flavor of the arguments, we enclose a lemma that 
is used many times in the proofs of the results above. 

Lemma. (Special curve lemma, [331 2.8]) Let E be a locally convex vector space. 
Let Xn be a sequence which converges fast to x in E; i.e., for each k G N the 
sequence n^{xn — x) is bounded. Then the infinite polygon through the can be 
parameterized as a smooth curve c : R —> A such that c(^) = Xn and c(0) = x. 

2.5 Remark 

Convenient calculus (i.e., having properties 6 and 7) exists also for: 

• Real analytic mappings; see (32] or [331 Chapter II]. 

• Holomorphic mappings; see [40] or [331 Chapter II] (using the notion of 
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• Many classes of Denjoy Carleman ultradifferentiable functions, both of Beurl- 
ing type and of Roumieu-type, see [35l [36l 1^1^ . 

• With some adaptations, Lip^; see [5S]. 

• With more adaptations, even (7^’“ (the fc-th derivative is Holder-continuous 
with index a); see EH], EZI- 

The following result is very useful if one wants to apply convenient calculus to 
spaces which are not tied to its categorical origin, like the Schwartz spaces S, T), 
or Sobolev spaces; for its uses see [50] and [38] . 

Theorem 2.1. [2H1 theorem 4.1.19] Let c : R —>■ i? fee a curve in a convenient 
vector space E. Let V <Z E' he a subset of bounded linear functionals such that 
the bornology of E has a basis of a{E,V)-closed sets. Then the following are 
equivalent: 

1. c is smooth 

2. There exist locally bounded curves : R —> i? such that I o c is smooth 
R —>■ R with {£ o = i o c^, for each £ €V. 

If E is reflexive, then for any point separating subset V C E' the bornology of 
E has a basis of cr(E, V)-closed subsets, by [521 4.1.23]. 


3 Manifolds of mappings and regular Lie groups 

In this section I hope to demonstrate how convenient calculus allows for very short 
and transparent proofs of the core results in the theory of manifolds of smooth 
mappings. 


3.1 The manifold structure on C°°{M,N) 

Let M be a compact (for simplicity’s sake) finite dimensional manifold and N a 
manifold. We use an auxiliary Riemann metric g on N. Then 


zero se^ion^...^ Ow 

TN ^ 

open 


(TTAT.exp®) 



yNxNQ 


open 


N X N 


C°°{M,N), the space of smooth mappings M ^ N, has the following manifold 
structure. A chart, centered at / G C°°{M,N), is: 

C“(M, N)Z)Uf = {g-. (f,g){M) c Uf C r(/*T7V) 

Ufig) = (7r7v,expS)-i o (/,g), Uf{g){x) = {exp^^y^)-\g{x)) 
{uf)-\s) = exp^os, (M/)-i(s)(a::) = exp®^^^ (s(a:)) 

Lemma 3.1. ^“(R, r(M; f*TN)) = r(R x M; pra* f*TN) 
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This follows by cartesian closedness after trivializing the bundle f*TN. 
Lemma 3.2. The chart changes are smooth (C°^) 

3 s {ttn, exp®) o s 1 -^ (ttat, exp®)“^ o (/ 2 , exp^^ os) 

Since they map smooth curves to smooth curves. 

Lemma 3.3. C'°“(R, TV)) ^ C'°°(R x M,iV). 

By lemma [ST] 

Lemma 3.4. Composition C°°{P,M) x C°°{M^N) C°^{P,N), {f,g) ^ go f, 

is smooth 

Since it maps smooth curves to smooth curves. 

Corollary 3.5. The tangent bundle of C°°{M, N) is given by 

TC°°{M,N) = C°°{M,TN) > C^{M,N). 

This follows from the chart structure. 

3.2 Regular Lie groups 

We consider a smooth Lie group G with Lie algebra q = T^G modelled on con¬ 
venient vector spaces. The notion of a regular Lie group is originally due to 
[Sllil ESI EnilSZlISH] for Frechet Lie groups, was weakened and made more 
transparent by m, and then carried over to convenient Lie groups in [34) . see 
also ESI 38.4]. We shall write /i : G x G —^ G for the mutiplication with 
x.y = pi{x,y) = pix{y) = for left and right translation. 

A Lie group G is called regular if the following holds: 

• For each smooth curve X G G°“(]R, g) there exists a curve g G G°°(M.,G) 
whose right logarithmic derivative is X, i.e., 

fg(0) =e 

\dtg{t) =Te{fia(t))Xit)=X{t).g{t) 

The curve g is uniquely determined by its initial value (?(0), if it exists. 

• Put evolQ(Al) = 5 ( 1 ) where g is the unique solution required above. Then 
evolg : G°“(M, g) —> G is required to be G°° also. We have Evol^ := g{t) = 
evolcitX). 

Up to now, every Lie group modeled on convenient vector spaces is regular. 

Theorem 3.6. For each compact manifold M, the diffeomorphism group DifF(M) 
is a regular Lie group. 
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Proof. Diff(M) M). Composition is smooth by restriction. In¬ 

version is smooth: If t i— f{t, ) is a smooth curve in Diff(M), then f{t, 
satisfies the implicit equation f{t, f{t, )~^{x)) = x, so by the finite dimensional 
implicit function theorem, (t,x) f{t, )~^{x) is smooth. So inversion maps 
smooth curves to smooth curves, and is smooth. 

Let X(t,x) be a time dependent vector field on M (in X(M))). Then 

(t, x) = (t + s, Evol^(f, x)) satisfies the ordinary differential equation 

dt Evol(t, x) = X{t, Evol(t, x)). 

If X{s,t,x) G (^“(IR^, X(M)) is a smooth curve of smooth curves in X(M), then 
obviously the solution of the equation depends smoothly also on the further vari¬ 
able s, thus evol maps smooth curves of time dependant vector fields to smooth 
curves of diffeomorphism. □ 

Groups of smooth diffeomorphisms on R" 

If we consider the group of all orientation preserving diffeomorphisms Diff(]R") of 
R", it is not an open subset of C'°°(R",R"') with the compact (7°°-topology. So it 
is not a smooth manifold in the usual sense, but we may consider it as a Lie group 
in the cartesian closed category of Erolicher spaces, see [231 Section 23], with the 
structure induced by the injection / i-^- (/, f~^) G C'°“(R",R") x C'°“(R",R"). Or 
one can use the setting of ‘manifolds’ based on smooth curves instead of charts, 
with lots of extra structure (tangent bundle, parallel transport, geodesic structure), 
described in [321 US]; this gives a category of smooth ‘manifolds’ where those 
which have Banach spaces as tangent fibes are exactly the usual smooth manifolds 
modeled on Banach spaces, which is cartesian closed: N) and Diff(M) are 

always ‘manifolds’ for ‘manifolds’ M and N, and the exponential law holds. 

We shall now describe regular Lie groups in Diff(R") which are given by dif¬ 
feomorphisms of the form / = Mr -\-g where g is in some specific convenient vector 
space of bounded functions in (^““(R",R"). Now we discuss these spaces on R", 
we describe the smooth curves in them, and we describe the corresponding groups. 
These results are from m and from [SSJ [ 32 ] for the more exotic groups. 

The group Diffg(R") 

The space S(R”) (called I?ioo(R") by [62]) consists of all smooth functions which 
have all derivatives (separately) bounded. It is a Frechet space. By [64], the space 
S(R") is linearly isomorphic to £°°^s for any completed tensor-product between 
the projective one and the injective one, where s is the nuclear Frechet space of 
rapidly decreasing real sequences. Thus B(R"') is not reflexive, not nuclear, not 
smoothly paracompact. 

The space (^“(R, B(R”)) of smooth curves in S(R") consists of all functions c G 
C'°“(R"“''^, R) satisfying the following property: 

• For all k G N>o, a G N>q and each t G R the expression d^dfc{t,x) is 
uniformly bounded in x G R", locally in t. 
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To see this use Theorem [Q for the set {ev^^ : cc € K.} of point evaluations in 
S(R"). Here 9“ = and c'=(t) = d^f{t, ). 

Diffg(IR") = {/ = Id+(/ : g £ S(]R")”,det(I„ + dg) > e > O} denotes the 
corresponding group, see below. 

The group Diffw^.p (R”) 

For 1 < p < oo, the space = nfe>i the intersection of all LP- 

Sobolev spaces, the space of all smooth functions such that each partial derivative 
is in LP. It is a reflexive Frechet space. It is called X>lp(R”) in [S2- By [53], 
the space IF°“’P(]R") is linearly isomorphic to £p^5. Thus it is not nuclear, not 
Schwartz, not Montel, and smoothly paracompact only if p is an even integer. 

The space C°^(M., (MT)) of smooth curves in IF“’^’(]R") eonsists of all functions 

c £ M) satisfying the following property: 

• For all k £ N>o, a £ N"g the expression \\dfd^f(t, )|| is locally 

bounded near each t £ R. 

The proof is literally the same as for S(R"'), noting that the point evaluations 
are continuous on each Sobolev space with k > ^. 

Diff]^oo.p(R") = {f = Id+g : g G IF°“’^(R”)”,det(I„ + dg) > 0} denotes the 
corresponding group. 

The group Diff 5 (R”) 

The algebra S{R^) of rapidly decreasing functions is a reflexive nuclear Frechet 
space. 

The space C°°(R, 5(R")) of smooth curves in 5(R") consists of all functions c £ 
(^oo(]^n+i, gdiigfyi^g i/jg following property: 

• For all k,m G N>o and a £ N”q, the expression (1 + \x\’^)'^d^df:c{t,x) is 
uniformly bounded in x G R", locally uniformly bounded m t £ R. 

DiffJ(R") = {/ = Id+(? : 5 £ 5(R”)",det(I„ + dg) > 0} is the corresponding 
group. 

The group Diffc(R") 

The algebra C^(R") of all smooth functions with compact support is a nuclear 
(LF)-space. 

The space C°°(R, (^^(R")) of smooth curves in C“(R"') consists of all functions 
f G C°°(R"“''^, R) satisfying the following property: 

• For each compact interval [a, b] in R there exists a compact subset K C R" 
such that f{t,x) = 0 for {t,x) G [a,b] x (R" \K). 

DiffcCR") = {f = ^d+g : g G C^(R")",det(I„ + dg) > O} is the correponding 
group. The case Diffc(R") is well-known since 1980. 
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Ideal properties of function spaces 

The function spaces discussed are boundedly mapped into each other as follows: 

^ 5(K")-^ ^ S(R") 

and each space is a bounded locally convex algebra and a bounded S(]R"')-module. 
Thus each space is an ideal in each larger space. 

Theorem 3.7 ([^ and [38]). The sets of diffeomorphisms Diffc(R"'), DifT 5 (R"), 
Diff//oo (K"), and DiffB(]R”) are all smooth regular Lie groups. We have the fol¬ 
lowing smooth injective group homomorphisms 

Difrc(K”)-^ Diffs(R")-^ Diffw~.p(R”)-^ DiffB(R") . 

Each group is a normal subgroup in any other in which it is contained, in particular 
in DiffB(]R"'). 

The proof of this theorem relies on repeated use of the Faa di Bruno formula for 
higher derivatives of composed functions. This offers difficulties on non-compact 
manifolds, where one would need a non-commutative Faa di Bruno formula for 
iterated covariant derivatives. In the paper [38] many more similar groups are 
discussed, modeled on spaces of Denjoy-Carleman ultradifferentiable functions. It 
is also showm that for p > 1 the group Diff^v^.pnii (K”) is only a topological group 
with smooth right translations — a property which is similar to the one of finite 
order Sobolev groups Diff^y^.p (R"). Some of these groups were used extensively 
in |52) . 

Corollary 3.8. DifFB(R”) acts on Fc, Fg andVH°° of any tensorbundle over R” 
by pullback. The infinitesimal action of the Lie algebra Tb(R"') on these spaces by 
the Lie derivative thus maps each of these spaces into itself. A fortiori, Diff//oo (R") 
acts on Fs of any tensor bundle by pullback. 

3.3 Manifolds of immersions and shape spaces 

For finite dimensional manifolds M, N with M compact, Emb(M, iV), the space 
of embeddings of M into N, is open in C°°{M,N), so it is a smooth manifold. 
Diff(M) acts freely and smoothly from the right on Emb(M, TV). 

Theorem 3.9. Emb(M, A^) —Emb(M, TV)/Diff(M) = B{M,N) is a principal 
fiber bundle with structure group Diff(M). 

Note that B{M, N) is the smooth manifold of all submanifolds of N which are 
of diffeomorphism type M. Therefore it is also called the nonlinear Grassmannian 
in 130], where this theorem is extended to the case when M has boundary. From 
another point of view, B{M,N) is called the differentiable Chow variety in |43j . 
It is an example of a shape space. 
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Proof. We use an auxiliary Riemannian metric g on N. Given / G Emb(-M, N), we 
view f{M) as a submanifold of N and we split the the tangent bundle of N along 
f{M) as TN\f(^M) = Nor(/(M)) ®Tf[M). The exponential mapping describes a 
tubular neighborhood of /(M) via 

Nor(/(M)) Wf(M) f{M). 

li g ■. M ^ N is C^-near to /, then (.p{g) := f~^ opf{M) G DifF(M) and we may 
consider g o p{g)~^ € r(/*W/(M)) C r(/* Nor(/(M))). This is the required local 
splitting. □ 

Imm(M, iV), the space of immersions M ^ N, is open in C°°{M,N), and is 
thus a smooth manifold. The regular Lie group Diff(M) acts smoothly from the 
right, but no longer freely. 

Theorem 3.10 f [17)1. For an immersion f : M ^ N, the isotropy group 
Diff (M) f = {p& Diff (M) : / o (^ = /} 

is always a finite group, aeting freely on M; so M —^ M/ DifF(M) f is a covering 
of manifolds and f factors to f = f op. 

Thus N) —>■ Imm(M, TV)/DifF(M) =: Bi{M,N) is a projection onto 

an infinite dimensional orbifold. 

The space Bi{M,N) is another example of a shape space. It appeared in the 
form of Bi{S^,'Sf), the shape space of plane immersed curves, in [35] and [15] . 


4 Weak Riemannian manifolds 

If an infinite dimensional manifold is not modeled on a Hilbert space, then a 
Riemannian metric cannot describe the topology on each tangent space. We have 
to deal with more complicated situations. 

4.1 Manifolds, vector fields, differential forms 

Let M be a smooth manifold modeled on convenient vector spaces. Tangent vectors 
to M are kinematic ones. 

The reason for this is that eventually we want flows of vector fields, and that 
there are too many derivations in infinite dimensions, even on a Hilbert space FI: 
Let a G F(fFl, H)he a, continuous linear functional which vanishes on the subspace 
of compact operators, thus also on H ®H. Then / i->- a{d‘^f{0)) is a derivation at 
0 on since a{d^{f .g){Q)) = a^d"^ f (0).g{0) + df (0) dg{0) + dg{0) df (0) + 

f{0).d^g{0)) and a vanishes on the two middle terms. There are even non-zero 
derivations which differentiate 3 times, see |33l 28.4]. 

The (kinematic) tangent bundle TM is then a smooth vector bundle as usual. 
Differential forms of degree k are then smooth sections of the bundle Lgj(.g^(TM; K) 
of skew symmetric fc-linear functionals on the tangent bundle, since this is the 
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only version which admits exterior derivative, Lie derivatives along vector held, 
and pullbacks along arbitray smooth mappings; see [33l 33.21]. The de Rham 
cohomology equals singular cohomology with real coefficients if the manifold is 
smoothly paracompact; see [531 Section 34]. If a vector held admits a how, then 
each integral curve is uniquely given as a how line; see [331 32.14]. 

4.2 Weak Riemannian manifolds 

Let (M, g) be a smooth manifold modeled on convenient locally convex vector 
spaces. A smooth Riemannian metric on M is called weak if ■ T^M T*M 
is only injective for each x S M. The image g{TM) C T*M is called the smooth 
cotangent bundle associated to g. Then g~^ is the metric on the smooth cotan¬ 
gent bundle as well as the morphism g{TM) —>■ TM. We have a special class 
of 1-forms flg{M) := T{g{TM)) for which the musical mappings makes sense: 
a'^ = g~^a € X(M) and = gX. These 1-forms separate points on TM. The 
exterior derivative is dehned by d : M^g{M) —>• Vi^{M) = since 

the embedding g(TM) C T*M is a smooth hber linear mapping. 

Existence of the Levi-Civita covariant derivative is equivalent to: The metric 
itself admits symmetric gradients with respect to itself. Locally this means: If M 
is c°°-open in a convenient vector space Vm- Then: 

D^^x 9 x{X,Y) = g^{X,grad^g{x){X,Y)) = g^igrad^ g{x){X, X),Y) 

where denote the directional derivative at x in the direction X, and where 

the mappings gradj^ g and symgrad 2 g : M x Vm x Vm —>■ Vm, given by (x, X) i—>- 
gradi 2 g{x){X, X), are smooth and quadratic in A € Vm- 

4.3 Weak Riemannian metrics on spaces of immersions 

For a compact manifold M and a finite dimensional Riemannian manifold {N, g) we 
can consider the following weak Riemannian metrics on the manifold lmm(M, N) 
of smooth immersions M ^ N: 

G°f{h,k)= [ g{h,k)vo\{f*g) 

Jm 

G){h,k)= [ g((l + A/*5)^/i,fc)vol(rg) 

Jm 

G^h,g)= [ $(/)g(/i,fc)vol(rg) 

Jm 

Here vol{f*g) is the volume density on M of the pullback metric g = f*g, where 
is the (Bochner) Laplacion with respect to g and g acting on sections of 
f*TN, and where $(/) is a positive function of the total volume Vol{f*g) = 
Im vol{f*g), of the scalar curvature Scal(/*g), and of the mean curvature Tr(5''^), 
being the second fundamental form. See [3] , [TU] for more information. All these 
metrics are invariant for the right action of the reparameterization group Diff(M), 
so they descend to metrics on shape space Bi{M,N) (off the singularities) such 


the L^-metric, 
a Sobolev metric of order s, 
an almost local metric. 


Manifolds of mappings and shapes 


13 


that the projection Imm(M, A^) —is a Riemannian submersion of a 
benign type: the G-orthogonal component to the tangent space to the DifF(M)- 
orbit consists always of smooth vector fields. So there is no need to use the notion 
of robust weak Riemannian metrics discussed below. 

Theorem 4.1. 1. Geodesic distance on Imm(M, iV), defined as the infimum of 
path-lenghts of smooth isotopies between two immersions, vanishes for the Lf- 
metric G°. 

2. Geodesic distance is positive on Bi{M,N) for the almost local metric G^ if 
‘h(/) > 1 + ATr(S''^), or if $(/) > AYo\{f*g), for some A> 0. 

3. Geodesic distance is positive on Bi{M, N) for the Sobolev metric G® i/s > 1. 
f. The geodesic equation is locally well-posed on Imm(M, TV) for the Sobolev 

metric G® if s > 1, and globally well-posed (and thus geodesically complete) on 
Imm(S'^, M"), if s>2. 

1 is due to [IH] for Bi{S^,Mf), to for Bi{M,N) and for Diff(M), which 
combines to the result for Imm(M, A^) as noted in [3]. 2 is proved in [T^. For 3 
see [3]. 4 is due to m and m- 

4.4 Weak right invariant Riemannian metrics on regular Lie 
groups 

Let G be a regular convenient Lie group, with Lie algebra g. Let /i : G x G —)• G 
be the group multiplication, the left translation and the right translation, 
Sx{y) = = xy = g{x,y). Let L,R : q ^ 3C{G) be the left- and right- 

invariant vector field mappings, given by Lx{g) = Te{yig).X and Rx = Te{pL^).X, 
respectively. They are related by Lx{g) = RAd{g)x{9)- Their flows are given by 

Flf-(5)=5.exp(tX) = ^-P(‘^)(5), 

Flf'^(5) = exp(tX).g = Mexp(tJf)(5)- 

The right Maurer-Cartan form k = k'^ £ n^(G,g) is given by Kx{ 0 ■= Tx{tG 
The left Maurer-Cartan form h) S n^(G,g) is given by Kx{C) '■= Tx{g,x-i) ■ C- 
satisfies the left Maurer-Gartan equation dn — = 0, where [ , ]a 

denotes the wedge product of g-valued forms on G induced by the Lie bracket. 
Note that k]a(C, ??) = 

k) satisfies the right Maurer-Gartan equation dn -\- /v]g = 0. 

Namely, evaluate dM on right invariant vector fields Rx, Ry for X,Y € g. 

{dK^){Rx,RY) = RxiK'-iRy)) - RviK'-iRx)) - k^{[Rx,Ry]) 

= Rx{Y) - Ry{X) + [X, F] = 0 - 0 + K(i?x), k''{Ry)]- 

The (exterior) derivative of the function Ad : G —>• GL(g) can be expressed by 

dAd = Ad.(ad OK*) = (ad ok’’). Ad 

since we have dAd(T^g.X) = ^|o Ad(( 7 . exp(tA)) = Ad{g). ad^K^Tpg.X)). 
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Let 7 = 0 xg—>-]Rbea positive-definite bounded (weak) inner product. Then 

').r]) =j{Ki^), K{ri)) 

is a right-invariant (weak) Riemannian metric on G and any (weak) right-invariant 
bounded Riemannian metric is of this form, for suitable 7 . Denote by 7 : g —>■ g* 
the mapping induced by 7 , from the Lie algebra into its dual (of bounded linear 
functionals) and by {a,X)g the duality evaluation between a G g* and X G q. 

Let g ■. \a,h] ^ G he a. smooth curve. The velocity field of g, viewed in the 
right trivializations, coincides with the right logarithmic derivative 

^’'( 5 ) = ^).dtg = K{dtg) = ig*K){dt). 

The energy of the curve g(t) is given by 

Eig) = ^j Jg{g',g')dt =^J -/{ig*K){dt), {g*K){dt)) dt. 

For a variation g{s,t) with fixed endpoints we then use that 

d{g*K){dt,ds) = dt{g*K{ds)) - ds{g*K{dt)) - 0, 


partial integration, and the left Maurer-Cartan equation to obtain 
dsE{g) = ^ j 2-f{dsig*K){dt), ig*K)idt)) dt 

= [ lidtig*K){ds) - d{g*K){dt,ds), {g*K){dt)) dt 


= - ( dt{g*K){dt)) dt 

J a 


- / l{[{9*E){dt),{g*K){ds)\, {g*K){dt))dt 


= - [ {l{dti9*K){dt)), {g*K){ds))^dt 

J a 


- / {i{{9*E){dt)), ad(g,^)(a^){g*K){ds))^dt 

J a 

= - [ {i{dtig*K){dt)) + {ad(^g-n)(dt)Tlii9*K){dt)), ig*K){ds))^ dt. 

J a 

Thus the curve g{0,t) is critical for the energy if and only if 

^idt{g*K)idt)) -f (ad(g.«)(a^))*7((g*«;)(at)) = 0. 

In terms of the right logarithmic derivative u : [a, 6 ] —s- g of g : [a, b] —> G, given 
by u{t) := g*K{dt) = Tg(t){p,^d) ) . g'{t), the geodesic equation has the expression 


dtu = —j ^ ad('u)* 7 (u) 
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Thus the geodesic equation exists in general if and only if ad(X)* 7 (X) is in the 
image of 7 : g —g*, i.e., 

ad(X)* 7 (X) e 7 (g) 

for every X € X. This condition then leads to the existence of the Christoffel 
symbols. Arnold in [5] asked for the more restrictive condition ad(X)* 7 (F) € 7 (g) 
for all X, T € g. The geodesic equation for the momentum p := 7 (m) is 

Pt = -ad(7"\p))V 

There are situations, see theorem l4.5l or [^, where only the more general condition 
is satisfied, but where the usual transpose ad^(X) of ad(X), 

ad^(X) := 7 “^ o ad^ o 7 


does not exist for all X. 

We describe now the covariant derivative and the curvature. The right trivi- 
alization (ttgjK’') : TG —>■ G x g induces the isomorphism R : G°°(G, g) —> X(G), 
given by R{X){x) := Rx{x) := Te{p^) ■ X{x), for X G G°°(G,g) and x G G. Here 
X(G) := r(TG) denotes the Lie algebra of all vector fields. For the Lie bracket 
and the Riemannian metric we have 

[i?x, = i?(-[A, y]0 + dY-Rx-dX-Rr), 

r-^[Rx,Ry] = -[X,y]g + Rx(r) - Ry(x), 

'yco(Rx(x),RY(x)) = j(X(x),Y(x)), x G G. 

In what follows, we shall perform all computations in C°°(G,g) instead of X(G). 
In particular, we shall use the convention 

VxY :^R-^(Vji^RY) for X, T G G“(G, g) 

to express the Levi-Civita covariant derivative. 

Lemma 4.2. 0 3.3] Assume that for all ^ G g the element ad(^)* 7 (^) G g* is in 
the image 0/7 : g —)• g* and that £, 1 —>■ 7 ~^ ad(^)* 7 (^) is bounded quadratic (or, 
equivalently, smooth). Then the Levi-Civita covariant derivative of the metric 7 
exists and is given for any X,Y G G°°(G, g) in terms of the isomorphism R by 

VxY = dY.Rx + p{X)Y - i ad(X)r, 

where 

PiOv = (adj^^ 7(^ + b) - ad|_^ 7(^ - 7)) = 57“^ (adj 7(77) + ad; 7(^)) 

is the polarized version, p : g —> L(g,g) is bounded, and we have p{£,)r] = p{r])£,. 
We also have 

= ^7(C,ad(p)C) + 

l{p{C)rh C) + l{p(.v)C, C) + 7(p(C)C, 0 = 0- 
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For X,Y e C“(G, g) we have 

[i?x, ad(y)] = ad(i?x(P)) and [Rx,piY)] = p{Rx{Y)). 

The Riemannian curvature is then computed by 

= [Vjf, Vy] - '\/_[x,Y]^+RxiY)-RY(X) 

= [Rx + Px — ^ adjc, Ry + py — ^ ady] 

- R{-[X,Y], + RxiY) - Ry{X)) - p{-[X,Y]^ + Rx{Y)- Ry{X)) 

+ i ad(-[X,y]g + Rx{Y)- Ry{X)) 

= [px,Py] +P[x,f]s - -j[px,a-<^Y] + -[py,adx] - - ad[x,y]j ■ 

which is visibly a tensor field. 

For the numerator of the sectional curvature we obtain 

7(7^(x, r)X, Y) = -fipxPYX, Y) - -/{PYPXX, Y) + 7(P[x.v]^, Y) 

- ^jipx [Y, X], r) + ^j{[Y, pxX], Y) 

+ 0 - ]p{[x, pyXIy) - i7([[x,y],X],y) 

= ^{pxX, pyY) - WpxYf^ + ^ II [X, y] ||2 

-^Jix, [Y,[X,Y]]) + ^^iY,[X,[X,Y]]) 

= 7(pxX, pyY) - WpxYf^ + ^II[X, y]||2 

-j{pxY,[X,Y]])+j{Y,[X, [X,Y]]). 

If the adjoint ad(X)^ : 0 “>■ 0 exists, this is easily seen to coincide with Arnold’s 
original formula [5], 

7(7^(A, y)X, y) = - J II ad(X)Ty + ad(y)^A||2 + 7 (ad(A)^A, ad(y)^y) 

+ i7(ad(A)Ty - ad(y)^A, ad(A)y) + ^ || [A, Y%. 

4.5 Weak right invariant Riemannian metrics on diffeomor- 
phism groups 

Let A be a finite dimensional manifold. We consider the following regular Lie 
groups: Diff(A), the group of all diffeomorphisms of A if A is compact. Diffc(A), 
the group of diffeomorphisms with compact support, if A is not compact. If 
A = R", we also may consider one of the following: DifF 5 (K"), the group of all 
diffeomorphisms which fall rapidly to the identity. Diffu/oo,p(IR."), the group of all 
diffeomorphisms which are modelled on the space iy°“’P(R")"', the intersection of 
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all VF^’^’-Sobolev spaces of vector fields. The last type of groups works also for 
a Riemannian manifold of bounded geometry {N,g); see |21j for Sobolev spaces 
on them. In the following we write Diffyv(A^) for any of these groups. The Lie 
algebras are the spaces of vector fields, where A G with 

the negative of the usual bracket as Lie bracket. 

A right invariant weak inner product on Diff_ 4 (iV) is given by a smooth positive 
definite inner product 7 on the Lie algebra which is described by the 

operator L = 7 : X^(iV) — 5 > X^(A^)' and we shall denote its inverse hy K = L~^ : 
L{Xa{N)) —>• Xa{N). Under suitable conditions on L (like an elliptic coercive 
(pseudo) differential operator of high enough order) the operator K turns out to 
be the reproducing kernel of a Hilbert space of vector fields which is contained in 
the space of either (bounded with respect to g) or vector fields. See 
[551 Chapter 12], [33], and [52] for uses of the reproducing Hilbert space approach. 
The right invariant metric is then defined as in 14.41 where ( , )xa{n) th® duality: 


G!f{Xo^,Yo^)= Gfd(^, Y) = 7 (A, F) = (L(X), F) 3 e^(^). 


Mis-using the notation for L we will often also write 

Gfd(X,F)= [ g{LX,Y)vol{g). 

J N 

Examples of metrics are: 

G°d(X,r)= [ g{X,Y)volig) 

J N 

Gfd(X,r)= [ 5 ((l + A®)«A,r)vol(g) 

J N 

Gf^[X,Y)= f X'.Y'dx = - [ X"Ydx 

Jr. Jr 

The geodesic equation on Diff^(A). As explained in 14.41 the geodesic equation is 
given as follws: Let (p : [a, h] DifF^(A) be a smooth curve. In terms of its right 
logarithmic derivative u : [a, &] ^ X^(A), u{t) := (p*K{dt) = p^'{t) o , the 

geodesic equation is 

L{ut) = L{dtu) = —eLd{u)*L{u). 

The condition for the existence of the geodesic equation is as follows: 

A K{&d{X)*L{X)) 

is bounded quadratic Xyv(A) —)• Using Lie derivatives, the computation 

of ad^ is especially simple. Namely, for any section uj oi T*N ^ vol and vector 
fields ^,7] G X_a_{N), we have: 


the metric, 
a Sobolev metric of order s, 
where X,Y G X^(IR.). 


/ [ (w,(??)) = - / (£j(w),ry), 

JN JN JN 
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hence ad^(a;) = +£^(a;). Thus the Hamiltonian version of the geodesic equation 
on the smooth dual L{Xa{N)) C (T*A^ (g> vol) becomes 

dta = — ad^(„) a = —CK{a)Oi, 
or, keeping track of everything, 


dt^p = uo (fi, dta = —Cua u = K{a) = a**, a = L{u) = tt'’. 

Theorem 4.3. Geodesic distance vanishes on Diff_ 4 (A^) for any Sobolev metric of 
order s<^. If N = x C with C compact, then geodesic distance vanishes also 
for s = i. It also vanises for the if-metric on the Virasoro group M x Diff^(]R). 

Geodesic distance is positive on Diff^(A'^) for any Sobolev metric of order s > 1. 
//dim(iV) = 1 then geodesic distance is also positive for s > ^. 

This proved in and [3]. Note that low order Sobolev metrics have 

geodesic equations corresponding to well-known non-linear PDEs: On DifF(S'^) or 
Diff^(]R) the L^-geodesic equation is Burgers’ equation, on the Virasoro group it 
is the KdV equation, and the (standard) i/^-geodesic is (in both cases a variant 
of) the Camassa-Holm equation; see [71 7.2] for a more comprehensive overview. 
All these are completely integrable infinite dimensional Hamiltonian systems. 

Theorem 4.4. Let {N, g) be a compact Riemannian manifold. Then the geodesic 
equation is locally well-posed on Diff^(V) and the geodesic exponential mapping 
is a local diffeomorphism for a Sobolev metric of integer order s > 1. For a 
Sobolev metric of integer order s > geodesic equation is even globally 

well-posed, so that (DifF^(V), G®) is geodesically complete. This is also true for 
non-integer order s if N = K". 

For N = S^, the geodesic equation is locally wellposed even for s > |. 

For these results see 0, m, m, m- 

Theorem 4.5. [B] For A € {Cf°let 

Ai(R) = {/ G C°°{R) : /' G Am , /(-oo) = 0} 

and Zet Diff^j (R) = {(/? = Id-|-/ : / G Ai(M),/'> —1}. These are all regular Lie 
groups. The right invariant weak Riemannian metric G^ (X, Y) = X'Y' dx 
is positive definite both on Diff_ 4 (R) where it does not admit a geodesic equa¬ 
tion (a non-robust weak Riemannian manifold), and on DifF^j(R) where it admits 
a geodesic equation but not in the stronger sense of Arnold. On Diff_ 4 ,^ (K) the 
geodesic equation is the Hunter-Saxton equation 

I rx 

(V5t) o = n, Ut = -uuj; + - {uoo{z)f dz , 


and the induced geodesic distance is positive. We define the R-map by: 

R : Diff^, m A(K, ir>-2) c Am R), Rm = 2 - 1) • 
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The R-map is invertible with inverse 

1 r 

: M(M,K>_2) Diffxi(R), R~^{j){x) = X + - dx . 

J — OO 

The pull-back of the flat L'^-metric via R is the -metric on Diff_ 4 (R), i.e., 
R*{•, •)(dx) = ■ Thus the space (Diff_ 4 j (R), iJ^) is a flat space in the 

sense of Riemannian geometry. There are explicit formulas for geodesics, geodesic 
distance, and geodesic splines, even for more restrictive spaces Ai like Denjoy- 
Carleman ultradifferentiable function classes. There are also soliton-like solutions. 
(Diff^j (R), ) is geodesically convex, but not geodesically complete; the geodesic 

completion is the smooth semigroup Mon^j^ = {ip = Id+/ : / £ Mi(R), f > —!}• 
Any geodesic can hit the subgroup Diff^(R) C Diff^j(R) at most twice. 

5 Robust weak Riemannian manifolds and Rie¬ 
mannian submersions 

Another problem arises if we want to consider Riemannian submersions, in par¬ 
ticular shape spaces as orbits of diffeomorphism groups, as explained in m- 

5.1 Robust weak Riemannian manifolds 

Some constructions may lead to vector fields whose values do not lie in T^M, but 
in the Hilbert space completion T^M with respect to the weak inner product gx- 
We need that Uxgm forms a smooth vector bundle over M. In a coordinate 
chart on open U C M, TM\u is a trivial bundle U xV and all the inner products 
gx,x £ U define inner products on the same topological vector space V. They all 
should be bounded with respect to each other, so that the completion V of V with 
respect to gx does not depend on x and [_]x^u = U x V. This means that 

forms a smooth vector bundle over M with trivialisations the linear 
extensions of the trivialisations of the bundle TM ^ M. Chart changes should 
respect this. This is a compatibility property between the weak Riemannian metric 
and some smooth atlas of M. 

Definition A convenient weak Riemannian manifold (M, g) will be called a 
robust Riemannian manifold if 

• The Levi-Civita convariant derivative of the metric g exists: The symmetric 
gradients should exist and be smooth. 

• The completions TxM form a smooth vector bundle as above. 

Theorem 5.1. If a right invariant weak Riemannian metric on a regular Lie 
group admits the Levi-Civita covariant derivative, then it is already robust. 

Proof. By right invariance, each right translation TpA extends to an isometric 
isomorphims TxG —> TxgG. By the uniform boundedness theorem these isomor¬ 
phisms depend smoothly on g £ G. □ 
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5.2 Covariant curvature and O’Neill’s formula 

In [431 2.2] one finds the following formula for the numerator of sectional curvature, 
which is valid for closed smooth 1-forms a, fd G on a weak Riemannianian 

manifold (M, g). Recall that we view g : TM ^ T*M and so g~^ is the dual inner 
product on g{TM) and a® = g~^{a). 

g{R{a^, , fi'^) = 

- \a^a\\\P\\l-i) - i/3“/3#(||a||2_i) -p + P^J)g-^{a, j3) 

(last line = —a**/3([a**,/3**]) -|-/3**a([a**,/?**]])) 

+ !I|[«^/3“]C 

This is called Mario’s formula since Mario Micheli derived the coordinate version 
in his 2008 thesis. Each term depends only on g~^ with the exception of the last 
term. The role of the last term (which we call the O’Neill term) will become 
clear in the next result. Let p : {E,gE) —>■ {B^gs) be a Riemannian submersion 
between infinite dimensional robust Riemann manifolds; i.e., for each b G B and 
X G Ef, := p~^{b) the tangent mapping T^p : (T^E^gE) {TbB,gB) is a surjective 
metric quotient map so that 

liailsB := inf{||Xx|lsE : X,, G T,E,T,p.X, = Cb}. 

The infinimum need not be attained in T^E but will be in the completion T^E. 
The orthogonal subspace {Y^ : gE{Xx,Tx{Eb)) = 0} will therefore be taken in 
Tx{Eb) in T^E. If ab = ) G gsiTbE) C T^B is an element in the gs- 

smooth dual, then p*ab := {Txp)*{ab) = gBiol’T^p ) : T^E —^ K is in T*E but 
in general it is not an element in the smooth dual gEiT^E). It is, however, an 
element of the Hilbert space completion gEiT^E) of the (/^-smooth dual gE{TxE) 
with respect to the norm || Ilg-G and the element g~^^{p*oib) ='■ {p*ab)‘^ is in the 

II 11 ge -completion T^E of T^E. We can call g^^(p*ab) =: {p*ab)^ the horizontal 
lift of al = gB^{ab) G TbB. 

Theorem 5.2. [?21 2.6] Let p : {E,gE) {B^pb) be a Riemann submersion 
between infinite dimensional robust Riemann manifolds. Then for closed 1-forms 
a, P G Ltg^{B) O’Neill’s formula holds in the form: 

gB{R^{a'^,P'^)f3Ka'^) = 9 e{R^{{ p*a)^, {p*0){p*0, 

+ |l|[(p*«)“,(p*/3)“]niL 

Proof. The last (O’Neill) term is the difference between curvature on E and the 
pullback of the curvature on R. □ 

5.3 Semilocal versiou of Mario’s formula, force, aud stress 

In all interesting examples of orbits of diffeomorphisms groups through a template 
shape, Mario’s covariant curvature formula leads to complicated and impenetrable 
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formulas. Efforts to break this down to comprehensible pieces led to the concepts 
of symmetrized force and (shape-) stress explained below. Since acceleration sits in 
the second tangent bundle, one either needs a covariant derivative to map it down 
to the tangent bundle, or at least rudiments of local charts. In [33] we managed 
the local version. Interpretations in mechanics or elasticity theory are still lacking. 

Let (M, g) be a robust Riemannian manifold, x G M, a, [3 G gxiT^M). Assume 
we are given local smooth vector fields Xa and X,g such that: 

1. Xaix) = aHx), Xfj^x) = I3*{x), 

2. Then — Xa is zero at x hence has a well defined derivative Dx{a^ — Xa) 
lying in Hom(T^,M, Ta;M). For a vector field Y we have Dx{a^ — Xa)-Yx = 
[Y, a'^ — Xa]{x) = £y (a* — Xa)\x- The same holds for /3. 

3. Cx^{a) = Cx^iP) = Cxf,{a) = Cxf,{P) = 0, 

4. [Xa,Xp\=Q. 

Locally constant I-forms and vector fields will do. We then define: 

X(a,P) : = ^d{g~^{a, P)), a I-form on M called the force, 

V{a, P){x) : = Dx{P^ - Xf 3 ).a'^{x) 

= d{P'^ — e TxM called the stress. 

V{a,P){x) -V{P,a){x) = [a^,P'^]{x) 

Then in terms of force and stress the numerator of sectional curvature looks as 
follows: 

g[R{a\P^)P^,a^){x) = Rn + i?i 2 -f i ?2 + .^3 , where 

i?n = \ {C\^{g-^){P,P) -2Cx^Cx,{g-^){a,p) Y C\^{g-^){a,a)){x), 

Ri 2 = (J-(a, a),V{p, p)) + {T{p, p),V{a, a)) - {T{a, P),V{a, p) + V{P, a)), 
R2 = {\\Ria,P)\\l., - {Ria,a)),TiP,P))^_,)ix), 

Rs = -lma,P)-V{p,a)\\l. 

5.4 Landmark space as homogeneus space of solitons 

This subsection is based on 132; the method explained here has many applica¬ 
tions in computational anatomy and elswhere, under the name LDDMM (large 
diffeomorphic deformation metric matching). 

A landmark q = {qi,..., qx) is an TV-tuple of distinct points in R”; landmark 
space Land^(R") C (R”)^ is open. Let = {q^,... ,q^) be a fixed standard 
template landmark. Then we have the surjective mapping 

evgo : Diff^(R") ^ Land^(R”), 

(fi ^ eVqo{g}) = (p{q°) = (i^((??),.. .,(p{q%)). 
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Given a Sobolev metric of order s > ^ + 2 on Diff^(R"), we want to induce a 
Riemannian metric on Land'^(]R") such that evqg becomes a Riemannian submer¬ 
sion. 

The fiber of ev^o over a landmark q = (po(Q^) is 

{ip e Diffx(IR") : p{q°) = q] = o {p ^ Diff^(IR”) : p{q°) = g°} 

= {p & Diff^(R") : p{q) = q} o p^ . 


The tangent space to the fiber is 


{Xopo-.X& X 5 (K"),X(g,) = 0 for all i}. 


A tangent vector Y o pQ £ Diff 5 (]R”) is -perpendicular to the fiber over q 
if and only if 


/ {LY,X)dx 

JR" 


0 yX with X{q) = 0. 


If we require Y to be smooth then F = 0. So we assume that LY = P^.5q^, a 

distributional vector field with support in q. Here Pi € But then 


Y{x) = L ^f'^P,.Sq.) = [ K{x-y)'^P^.Sq,{y)dy = '^K{x-qi).Pi 

i ^ i i 

Tq,q{eVqo).{Y o p^) = Y{qk)k = '^{K{qk - qi).Pi)k ■ 


Now let us consider a tangent vector P = (Pk) G Tg Land'^(M”). Its horizontal 
lift with footpoint po is o pg where the vector field on M" is given as 
follows: Let K~^{q)ki be the inverse of the {N x A^)-matrix K{q)ij = K{qi — qj). 
Then 


P^^’^ix) = ^K{x - q,)K-\q),,Pq , 

L(ph”(a:)) = - q.)K-\q),,Pq . 

'id 


Note that P^” is a vector field of class 

The Riemannian metric on Land'^ induced by the (/^-metric on Diff 5 (]R") is 

g^iP, Q) = G^„(ph”, g>'”) = [ (L(ph”), g''”) dx 

= f (^6{x - qi)K~^{q),jPj,^K{x - qk)K~^{q)kiQi) dx 
i,j k,i ' 

= Yl, ^~^id)ijKiq^-qk)K-^{q)kl{Pj,Ql) 

g!^{P,Q)=YK-Hd)ki{Pk,Qi). 


kA 
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The geodesic equation in vector form is: 


= -\^ ^ ^{q)kigT:&dK{qi-qj){K{q)in-K{q)jn)K '^{q)ji{qk,qi) 


qn 




+ ^{q)ki(^gradK{qi -qn),qi - qn'jqk ■ 


k^i 

The cotangent bundle T*Land^(]R") = Land'^(]R”) x ((R”)^)* 9 (<?, a). We shall 
treat M" like scalars; ( , ) is always the standard inner product on R". 

The metric looks like 

=J2K{q)ij{ai,l3j), K{q),j =K{qi-qj). 

'id 

The energy function is 

E{q,a) = ^{g^)~'^{a,a) = ^'^K{q),j{a^,aj) 

'id 

and its Hamiltonian vector field (using E^-valued derivatives to save notation) is 

N 


HE{q,a)= y] (^K{qk - q^)ai-^ + gradK{qi - qk){ai,ak)-^y 


d 


i,k—l 


' dqk 


dak^ 


So the Hamitonian version of the geodesic equation is the flow of this vector field: 

iqk =Y.iE{qi- qk)ai 

\dfe =-X]*gradii:(gi - gfc)(ai,afc) 

We shall use stress and force to express the geodesic equation and curvature: 
d = '^Kiqk - qi)a^, = '^K{qk - qi){ai, -^) 

i i,k 

X>(a,/?) : = '^dK{qi - gj)(af - stress. 




T){a, j3) — T’(/3, a) = = [a**, /3**], Lie bracket. 

Ei{a,l3) = i y] gradKfe -%)({«*,/3fe) + (ft,afc)) 

k 

iF{a,l3) : = '^{Fi{a, (3),dqi) = tdg~^{a,(3) the force. 


The geodesic equation on T* Land^^ (M^) then becomes 

iq = 

Id = —iF{a, a). 
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Next we shall compute curvature via the cotangent bundle. From the semilocal 
version of Mario’s formula for the numerator of the sectional curvature for constant 
1-forms a,/3 on landmark space, where K{qk — qi)o'i, we get directly: 

- {Via, a), Tip, 13)) - {ViP, P),Tia,a)) 

- 5 X] /?“ - /3])(ai: aj) 

- 2d‘^Kiq, - qj)ipj - p],a\ - q:“)(A, a^) 

+ dpKiq, - qj)ia\ - a],a\ - a*)(/3,, Pj)'^ 

- \\Tia,P)\\l.^+g-\Tia,a),TiP,P)). 

5.5 Shape spaces of submanifolds as homogeneous spaces 
for the diffeomorphism group 

Let M he & compact manifold and (iV, g) a Riemannian manifold of bounded ge¬ 
ometry as in subsection 14.51 The diffeomorphism group Diffyi(iV) acts also from 
the left on the manifold of Emb(M, N) embeddings and also on the non-linear 
Grassmannian or differentiable Chow variety BiM,N) = Emb(M, iV)/Diff(M). 
Eor a Sobolev metric of order s > -|- 2 one can then again induce a Rieman¬ 

nian metric on each Diff_ 4 (iV)-orbit, as we did above for landmark spaces. This 
is done in m, where the geodesic equation is computed and where curvature is 
described in terms of stress and force. 
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